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Abstract 

We investigate the occurrence of power terms in the running QCD 
coupling Ois{p) by analysing non-perturbative measurements at low 
momenta {p > 2GeV) obtained from the lattice three-gluon vertex. 

Our exploratory study provides some evidence for power contribu¬ 
tions to as{p) proportional to 1/p^. Possible implications for physical 
observables are discussed. 


1 Introduction 

The standard procedure to parametrise non-perturbative QCD effects in 
terms of power corrections to perturbative results is based on the Opera¬ 
tor Product Expansion (OPE). In this framework, the powers involved in 
the expansion are expected to be uniquely fixed by the symmetries and the 
dimension of the relevant operator product. It should be noted that, due to 
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the asymptotic nature of QCD perturbative expansions, power corrections 
are reshuffled between operators and coefficient functions in the OPE [jl|. 

The above picture has recently been challenged It was pointed 

out that power corrections which are not a priori expected from the OPE 
may in fact appear in the expansion of physical observables. Such terms 
may arise from (UV-subleading) power corrections to as{p), corresponding 
to non-analytical contributions to the /3-function. To illustrate this point, 
consider for example a typical contribution to a condensate of dimension 2a: 


dp^ 

pA2 


Q2) “-(^2) • 


( 1 ) 


2 A 2 

A power contribution to of the form would generate (from 

the UV limit of integration) a contribution to the condensate proportional to 
{^Y■ The fact that the dimension of such a term would be independent of a 
indicates that this contribution would be missed in a standard OPE analysis. 

Note that in the above manipulations z could be in principle any (real) 
number. The value z = 1 may in fact play a special role (see the discussion in 
Section as it would result in ^ contributions to physical processes whose 
existence has been conjectured fcr a long time, mainly in the framework of 
the UV renormalon Q. 

Clearly, the existence of a-independent power corrections, if demon¬ 
strated, would have a major impact on our understanding of non-perturbative 
QCD effects and may affect QCD predictions for several processes. For ex¬ 
ample, ^ contributions may be relevant for the analysis of r decays P, ||. 

Although the size of such corrections could in principle be estimated 
directly from experimental data, it would be highly desirable to develop a 
theoretical framework where the occurrence of these effects is demonstrated 
and estimates are obtained from first principles QCD calculations. Some 
steps in this direction were performed in []^, ^, where some evidence for an 
unexpected ^ contribution to the gluon condensate was obtained through 
lattice calculations. 


The aim of the present work is to test a method to detect the presence 
of power corrections in the running QCD coupling. Non-perturbative lattice 
estimates of the coupling at low momenta are compared with perturbative 
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formulae. Although at this stage our work is exploratory in nature and fur¬ 
ther simulations will be required to obtain a conclusive answer, our analysis 
provides some preliminary evidence for power corrections. The final goal is 
to investigate the possible link between OPE-independent power corrections 
to physical observables and power terms in the running coupling. 

The paper is organised as follows: in Section we briefly review some 
theoretical arguments in support of power corrections to as{p), illustrating 
the special role that may be played by ^ terms. In Section we explain 
the meaning of the lattice data and our strategy for the analysis. Some 
preliminary evidence for power corrections is discussed. Finally, in Section ^ 
we draw our conclusions. The appendix contains some technical details. 

2 Clues for ^ Corrections to as{p) 

Power corrections to as{p) can be shown to arise naturally in many physical 
schemes p|, The occurrence of such corrections cannot be excluded a priori 
in any renormalisation scheme. Clearly, given the non-analytic dependence 
of terms on power corrections cannot be generated or analysed 

in perturbation theory. In particular, the non-perturbative nature of such 
effects makes it very hard to assess their dependence on the renormalisation 
scheme, which is only very weakly constrained by the general properties of 
the theory. 

As discussed in the following, despite the arbitrariness a priori of the 
exponent z, several arguments have been put forward in the past to suggest 
that a likely candidate for a power correction to as{p) would be a term of 
order A^/p^, i.e. z = 1. 

2.1 Static Quark Potential and Confinement 

Consider the interaction of two heavy quarks in the static limit (for a more 
detailed discussion see [|^). In the one-gluon-exchange approximation, the 
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static potential V (r) can be written as 


V(r) (X as d'^k 


exp 


ik-r 




( 2 ) 


Clearly the above formula yields the Coulomb potential V{r) ^ 1/r. 
Using standard arguments of renormalon analysis, one may consider a gen¬ 
eralisation of @ obtained by replacing with a running coupling: 


V{r) oc / d'^kas{\k\^) 


2 , exp 


ik-r 


\k[- 


( 3 ) 


The presence of a power correction term of the form oc A^//c^ 

would generate a linear conhning potential V (r) ~ Kr. Note that a standard 
renormalon analysis of (|^) (see m for the details) reveals contributions to 
the potential containing various powers of r, but a linear contribution is 
missing. This is a typical result of renormalon analysis: renormalons can 
miss important pieces of non-perturbative information. 


2.2 An Estimate from the Lattice 


The lattice community has been made aware for some time of the possi¬ 
bility of non-perturbative contributions to the running coupling; for a clear 
discussion see . Consider the “force” dehnition of the running coupling: 


(^qq{Q) 


3,jdV{r} 
4 dr 



( 4 ) 


where again V{r) represents the static interquark potential. By keeping into 
account the string tension contribution to U(r), which can be measured in 
lattice simulations, one obtains a 1/Q^ contribution, whose order of mag¬ 
nitude is given by the string tension itself. Ironically, this term has been 
mainly considered as a sort of ambiguity, resulting in an indetermination in 
the value of a{Q) at a given scale. From a different point of view, such a 
term could be interpreted as a clue for the existence of a ^ contribution, 
and it also provides an estimate for the expected order of magnitude of it, 
at least in one (physically sound) scheme. 
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2.3 Landau pole and analyticity. 


It is well known that perturbative QCD formulae for the running of ag in¬ 
evitably contain singularities, which are often referred to as the Landau pole. 
The details of the analytical structure depend on the order at which the (3- 
function is truncated and on the particular solution chosen. The existence of 
an interplay between the analytical structure of the perturbative solution and 
the structure of non-perturbative effects has been advocated for a long time 
H- To illustrate this idea, consider the one-loop formula for the running 


coupling Q!s(p): 




bo log(^) 


(5) 


Here the singularity is a simple pole, which can be removed if one redehnes 
as{p) according to the following prescription: 

1 A2 


asip ) = 




+ 


6o(A2 — p'^ 
A2 


( 6 ) 


where a power correction of the asymptotic form p- appears. However, the 
sign of such a correction is the opposite of what one would expect from the 
results of [Q and from the considerations in Section p.l|, so that in the end 
one could envisage a more general formula for the regularised coupling: 




+ 


A2 


bo log(^) &o(A 2 -p 2 ) 


+ c 


p2 


(7) 


The message from (0) is that the perturbative coupling is not defined at the 
^ level, so the coefficient of the power correction is unconstrained, even after 
imposing the cancellation of the pole. 


At higher perturbative orders one encounters multiple singularities, which 
include an unphysical cut. There are several ways to regularise them. In 
particular, the method discussed in |T^ combines a spectral-representation 
approach with the Renormalization Group. The method was originally for¬ 
mulated for QED, but it has recently been extended to the QCD case |T^ . 


Other approaches can be conceived to achieve a systematic regularisation 
of the singularities arising from the Landau pole, order by order in perturba¬ 
tion theory. In this way one obtains formulae for as{p) that are well-dehned 
at all momentum scales. 
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Such formulae would be quite useful in the framework of our study, since 
power corrections are expected to be sizeable at scales close to the location of 
the Landau pole. However, for the purpose of the preliminary investigation 
discussed in the present paper, we shall limit ourselves to a simpler approach, 
where one tries to £t the data by simply adding power corrections to the 
perturbative expressions, without attempting a regularisation of the Landau 
pole. 


3 Lattice Data and Power Corrections 


3.1 




on the Lattice 


Several methods for computing as{p) non-perturbatively on the lattice have 
been proposed in recent years [^, 0,0, irj 0. In most cases, the goal of 
such studies is to obtain an accurate prediction for as{Mz), i.e. the running 
coupling at the Z peak, which is a fundamental parameter in the standard 
model. For this reason, lattice parameters are usually tuned as to allow the 
computation of as{p) at momentum scales of at least a few GeVs, where the 
two-loop asymptotic behaviour is expected to dominate and power contribu¬ 
tions are suppressed. However, the same methods can in principle be applied 
to the study of as{p) at lower momentum scales, where power-like terms may 
be sizeable. For this purpose, the best method is one where one can measure 
as{p) in a wide range of momenta from a single Monte Carlo data set. 

One method which fulhlls the above criterion is the determination of the 
coupling from the renormalised lattice three-gluon vertex function 


18, IS 


This is achieved by evaluating two- and three-point off-shell Green’s functions 
of the gluon held in the Landau gauge, and imposing non-perturbative renor- 
malisation conditions on them, for different values of the external momenta. 
By varying the renormalisation scale p, one can determine as{p) for differ¬ 
ent momenta from a single simulation. Obviously the renormalisation scale 
must be chosen in a range of lattice momenta such that both hnite volume 
effects and discretisation errors are under control. Such a dehnition of the 
coupling corresponds to a momentum-subtraction renormalisation scheme in 
continuum QCD [^. It should be noted that in this scheme the coupling is 
a gauge-dependent quantity. One consequence of this fact is that 1/p^ cor- 
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rections should be expected, based on OPE considerations. We will return 
to this issue when drawing our conclusions. 


The numerical results for as{p) that we use for our investigation were 
obtained from 150 conhgurations on a 16“^ lattice at /3 = 6.0. 

m 


For full details of the method we refer the reader to Ref. 


where 


such results were hrst presented. In order to detect violations of rotational 
invariance, different combinations of lattice vectors have sometimes been used 
for a hxed value of p^. This accounts for the graphical “splitting” of some 
data points. 


3.2 Models for Power Corrections 


As mentioned at the end of Section 2.3, in the present work we shall not ad¬ 
dress the general problem of dehning a regular coupling at all scales. For 
the purpose of our preliminary investigation, we shall compare the non- 
perturbative data for with simple models obtained by adding a power 
correction term to the perturbative formula at a given order. In order to 
identify momentum intervals where our ansatz hts the data, one should keep 
in mind that the momentum range should start well above the location of 
the perturbative Landau pole, but it should nonetheless include low scales 
where power corrections may still be sizeable. The requirement of keeping 
the effects of the hnite lattice spacing under control in the numerical data 
for as induces a natural UV cutoff on the momentum range. It is reassuring 
that intervals that fulhll these requirements can be identihed, as specihed in 
the following. 


One problem in this approach is the possible interplay between a de¬ 
scription in terms of (non-perturbative) power corrections and our ignorance 
about higher orders of perturbation theory. In particular, for the scheme that 
we consider, the three-loop coefficient of the /3-function is not known. Knowl¬ 
edge of such a coefficient would allow to perform a more reliable comparison 
of our estimates for the A parameter in our scheme with lattice determina¬ 
tions of A in a different scheme, for which the three-loop result is available 
In fact, although matching the A parameter between different schemes 


only requires a one-loop computation (because of asymptotic freedom), the 
reliability of such a comparison rests on the assumption that the value of A 
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in each scheme is fairly stable with respect to the inclusion of higher orders, 
which in turn implies that a sufficient number of perturbative orders has 
been considered in the dehnition of A. In practice, when working at two- 
or three-loop order, the value of A is still quite sensitive to the order of the 
calculation. For this reason, in the formulae for as{p) we shall append a 
subscript to the parameter A, to remind the reader that the value of such 
a parameter is expected to carry a sizeable dependence on the order of the 
perturbative calculation. 

Even within such limitations, in the following we will argue that it is 
possible to estimate the impact of three-loop effects in our model and that 
a description with power corrections seems to be stable with respect to the 
inclusion of such effects. 


3.3 Two-loop Analysis 

At the two-loop level, we consider the following formula: 


as{p) 


1 _ bi log(log(pVA|)) ^ 

bo log(pVA|) bo {bo log(pVA|))2 p 2 


( 8 ) 


By htting our data to (H) we obtain two sets of estimates for the pa¬ 
rameters {A 2 i,C 2 i), namely (0.84(1),0.31(3)) and (0.73(1),0.99(7)). The two 
results correspond to comparable values for Xdofy both cases we ob¬ 

tain Xdof — 1-8- In both cases, the momentum window extends up to p ~ 3 
GeV. We take the hrst set of values as our best estimate of the parameters 
as the corresponding value of A 21 is close to what is obtained from a “pure” 
two-loop £t, i.e. A 21 is stable with respect to the introduction of power cor¬ 
rections. Our choice for the value of A 21 will be a posteriori supported also 
by independent considerations at the three-loop level. The momentum range 
that we are able to describe (1.8 — 3.0 GeV) is fully consistent with what 
one would expect from general considerations based on the value of the UV 
lattice cut-off and the value of A 21 . Notice that choosing between the two 
sets of values makes quite a difference in the UV region, where power effects 
are largely suppressed. 

In summary, a two-loop description with power corrections based on (|^) 
hts well the data in a consistent momentum range. Our best £t of the data 
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to (U) is shown in Figure 1. We were also able to check that if one tries to 
determine the exponent z of the power correction from the £t, the best 

description of the data is obtained for z ^ 1. We interpret this result as a 
conhrmation of our theoretical prejudice z = 1. However, one should note 
that since that the quality of our data makes a full three-parameter £t very 
hard, the above check of the value of z and any other three-parameter £t that 
we mention in the following sections were in fact obtained by performing a 
very large number of two-parameter hts, corresponding to different (hxed) 
values of the third parameter. 



Figure 1: The best £t to (|). The crossed-circled points indicate the htting 
range. 
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3.4 Three-loop Analysis 


As already mentioned, a major obstacle for a three-loop analysisis is the fact 
that the hrst non-universal coefficient 62 of the perturbative /3-function is not 
known for our scheme. 

In order to gain insight, we start by performing a two-parameter £t to 
the standard three-loop expression for as{p), where the fitting parameters 
are A 31 and the unknown coefficient 62 - We call the £t estimate for 62 , 
to emphasise that we expect the effective value to provide an order of 
magnitude estimate of the true (unknown) coefficient 62 - Our best estimate 
for A 31 and 62 is A^i = 0.72(1), = 1.3(1), with Xdof ~ (see the dashed 

curve in Fig. 2). The error quoted for the £t parameters should always be 
interpreted within the effective description provided by the relevant formula. 

The momentum range where we obtain the best description of the data is 
p ~ 2 — 3 GeV. Our result for A^i provides (via perturbative matching) an es¬ 
timate for Ajjg, in very good agreement with the estimate in |^, which was 
obtained from the computation of the A parameter in a completely different 
scheme. Although both estimates are affected by our ignorance of higher 
loop effects, and our estimate also depends on the extra parameter bl-^-^, the 
agreement between the two results appears remarkable. In order to investi¬ 
gate the reliability of b^^^ as an estimate of 62 , we discuss in the appendix 
an argument which appears to provide a lower bound for the value of 62 in 
our scheme, namely 62 ^ 0.3. Our value for is therefore consistent with 
such a bound. 

Having obtained comparable values for Xdof froni the two-loop analysis 
with power corrections and from the “pure” three-loop analysis, one may 
be led to consider our results as evidence against the existence of power 
corrections, since so far they simply appear to provide an effective description 
of three-loop effects. 

However, we will argue now that there is room for power corrections even 
at the three-loop level. To this aim, consider the following three-loop formula 
with a power correction: 


1 bi log(L) 
boL boiboiy 
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+ C3Z—f-, (9) 

pZ 

where L = log(p^/A 3 ;) and is again to be determined from a fit. 

Fitting the data to (^, we obtain As^ = 0.72(1), = 1.0(1) and 

C 3 = 0.41(2), with Xdof ~ 1-8) in n momentum range 1.8 — 3.0 GeV (see Fig. 
2). The above result was in practice obtained by performing a large number 
of two-parameter fits for b^^^ and C 3 , for fixed values of A^i. The range of 
trial values for A^i was suggested by the results of the “pure” three-loop fit. 

We note the following: 

1 . the value for the scale parameter A^i is fully consistent with the previous 
determination from the “pure” three-loop description; 

2 . the value for is also reasonably stable with respect to the previ¬ 
ous determination and it is also consistent with the approximate lower 
bound for 62 discussed in the appendix; 

3. by comparing results from hts to (H) and (||), it emerges that 

C 2 AI 1 = 0.22(2) GeV^ ~ csAh = 0.21(2) GeVl (10) 

This approximate equality gives us conhdence in the presence of power 
corrections, as it indicates that the power terms providing the best fit to 
(H) and (I) are numerically equal. In other words, there appears to be 
no interplay between the indetermination connected to the perturbative 
terms and the power correction term, within the precision of our data, 
thus suggesting that a genuine ^ correction is present in the data. 

Finally, the coefficient of the power correction is of the order of magnitude 
expected from the arguments in sections 2.1 and 2 . 2 , that is, it is comparable 
to the standard estimate for the string tension squared. 

One may argue at this point that at the two-loop level we had to choose 
between two sets of values for (A 2 z,C 2 z), and that our choice is crucial for the 
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validity of (^). An a posterzon justification for our choice can be obtained 
from the following test: we plot a few values for as{p) as generated by the 
“pure” three-loop formula for A 3 ; = 0.72 and 62 = 1-0. Then, by htting such 
points to the “pure” two-loop formula, one gets A 21 ~ 0.84, i.e. the value for 
which ([l 0 |) holds. 


4 Conclusions 

We have discussed an exploratory investigation of power corrections in the 
running QCD coupling as{p) by comparing non-perturbative lattice results 
with theoretical models. Some evidence was found for 1 /p^ corrections, whose 
size was consistent with what is suggested by simple arguments from the 
static potential. 

At the technical level, our results need further conhrmation from the 
analysis of a larger data set and a study of the dependence of the £t param¬ 
eters on the ultraviolet and infrared lattice cutoff. Assuming our hndings 
are conhrmed at the technical level, one needs to address the issue of as¬ 
sessing the scheme dependendence of our results. As already discussed, the 
non-perturbative nature of power corrections makes it very hard to formu¬ 
late any theoretical procedure to estimate the impact of scheme dependence. 
The best one can do at this stage is to consider different renormalisation 
schemes and dehnitions of the coupling and gather numerical evidence and 
formal arguments supporting power corrections to as{p)- In this way, scheme- 
independent features may eventually be identified. For example, on the basis 
of our results, we note the following: 

• Theoretical arguments suggest 1/p^ corrections both for the coupling 
as dehned from the static potential and for the one obtained from the 
three-gluon vertex. The arguments for the former case were outlined in 
Sections 2.1 and 2 . 2 . As far as the coupling from the three-gluon vertex 
is concerned, 1/p^ corrections appear in an OPE analysis if one keeps 
into account the fact that such a coupling is a priori gauge dependent, 
so that a dimension 2 condensate appears in the relevant OPE. 

• In the static potential case, the theoretical argument also provides an 
estimate for the order of magnitude of the coefficient of the 1 /p^ correc- 
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tion, while in the three-gluon vertex case the OPE argument provides 
no estimate for it, suggesting instead that it may depend on the gauge. 
However, our numerical result in the Landau gauge is in striking agree¬ 
ment with the estimate for the static potential case. Although such an 
agreement may of course be accidental, it calls for further investiga¬ 
tion, which may be performed by attempting a similar calculation in a 
different gauge. 

The issue of scheme dependence will be the focus of our future work. 
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Appendix A 


Consider the perturbative matching between our scheme and the MS scheme 


a. 


— cums + CiOjjg -|- 02(62) -|- O(a^) 


As it is well known, c\ determines the ratio of the A parameters in the 
different schemes, while C 2 depends on ci and the difference between the 
value of 62 in our scheme and . We assume that at very high momentum 
values {p > 150 GeV) the running coupling follows the three-loop asymptotic 
formula. Then if one takes the value for A]^ from |^| and the value for 


^MOM scheme from the perturbative matching, the only unknown 

parameter in the above expression is the value of 62 in our scheme. By 
demanding that at the two-loop level the expansion of one coupling in powers 
of the other is still convergent (i.e. the convergence is better at two loops than 
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at one loop as the series are not yet displaying their asymptotic nature) we 
obtain an approximate lower bound for the unknown coefficient as 62 ^ 0.3. 
We have checked that such a technique provides sensible results for every 
couple of couplings for which a two-loop matching is known. 
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Figure 2: Fits to (|) (solid line) versus a pure three-loop fit (dashed line). 
The crossed-circled points are consistent with both ansatze, while the star- 
circled one is best htted by (^. 
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